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Abstract. Time-continuous non-anticipating quantum processes of nonde- 
molition measurements are introduced as the dynamical realizations of the 
causal quasi-measurements, which are described in this paper by the adapted 
operator-valued probability mcasures on the trajectory spaces of the gener- 
alized temporal observations in quantum open systems. In particular, the 
notion of physically realizable quantum filter is defined and the problem of its 
optimization to obtain the best a posteriori quantum state is considered. It 
is proved that the optimal filtering of a quantum Markovian Gaussian signal 
with the Gaussian white quantum noise is described as a coherent Markovian 
linear filter generalizing the classical Kalman filter. As an example, the prob- 
lem of optimal measurement of complex amplitude for a quantum Markovian 
open oscillator, loaded to a quantum wave communication line, is considered 
and solved. 



1. Introduction 

At present, due to invention of làser as a coherent quantum generator, the dy- 
namical problem of optimal reception of optical signals with quantum noise, taking 
into account the fundamental limitations caused by the quantum nature of the elec- 
tromagnètic waves, is becoming actual. Quantum information and communication 
theory which has been developing so far (See recent review yQ and the references 
in there) is based on static (nontemporal) , or single step (instantaneous) theory 
of quantum measurement, and quantum statistical inference does not take into 
account physical causality in due coarse of the dynamical propagation. As the 
result, the optimal estimators based on nontemporal quantum measurements are 
noncausal since they may depend at each intermediate time t not only on the past 
but also on future results of measurements. An attempt to develop the temporal, 
multistep variant of quantum measurement and decision theory to solve the prob- 
lems of the dynamical filtering in discrete time was made in [2J. However, because 
of the very restrictive class of quantum measurement process as a sequences of in- 
dependent, nonadaptive quantum measurements considered in this paper, even the 
simplest Gaussian problem for quantum linear filtering (quantum Kalman filtering) 
was not solved, with the exception of a degenerated case, and no time-continuous 
generalization is possible within this restricted method. 
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Here we show, while solving the filtering problems for a quantum diffusion, that 
these difíiculties can be avoided if a natural much wider class of temporal adaptive 
quantum measurements is considered, assuming their dependence not only on time, 
but also on the results of preceding measurements. Our dynamical model for the 
temporal quantum causal measurements has a natural time-continuous setup based 
on quantum stochastic differential equations. An appropriate quantum measure- 
ment device, that realizes the adaptive measurements, is called the nonanticipating 
(or causal) quantum filter with a memory. The optimization of estimation based 
on independent quantum measurements in the narrower class of filters without 
memory is in our treatment a problem with constraints leading to a poorer quality 
of filtration. As is proved below, the optimal time-continuous filtering of Markov 
Gaussian signals with the background quantum white noise is realized by a quantum 
linear filter based on the adaptive coherent quantum measurements. The optimal 
filter has a Markov memory such that the next optimal quantum measurement in 
general depends on the result of the present measurement but is independent of 
the preceding measurement results in this Markovian case. We show that such co- 
herent quantum filter can be realized by indirect heterodyne measurements adding 
an independent vacuum quantum noise, and by causal processing of the temporal 
measurement results using a classical linear Kalman-Bucy filter [3]. A nonlincar 
(and quasi-linear) generalization of this result for the non-Gaussian case, which 
corresponds to the classical nonlinear filtering theory of Stratonovich 4 , is also of 
our interest an will be considered elsewhere. 

Let us first consider a physical model of quantum open system with output 
channel, in which the problem of optimal quantum filtering arises naturally, giving 
a solution of this problem. Thcn we shall set up a general rigorous formulation of 
the filtering problem as the problem of quantum optimal temporal estimation, and 
derive the solution of this problem in the Gaussian case. 

2. Optimal observation of quantum oscillator at the output of 



Let us consider an electromagnètic oscillator with frequeney Í2 loaded on a pair 
of communication channcls, the transmission lines which are described by wave 
conductivity G and wave resistance R respectively, as shown in the figure, wherc 



O = 1/y/LC, G = y/ci/h, R — y/fa/c2, and propagation velocities of the waves 

Vi = í/y/lid, i =1,2 are assumed to be equal: v% — v 2 = v : 
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We consider the case when the lines are homogeneous and that the measuring 
device (receiver) which is set at their outputs of the lines (0 — 1) and (0 — 2) is 
ideally conjugated with these lines such that there is no reflection of the received 
information (i.e. incoming waves J_,V_) into the radiated noise (i.e. outgoing 
waves J+jV-if-) at the end of the lines. Such system is open, but it is entirely 
described by the dynamic variables of voltage V{t) and current J(t) on the contour 
(L, C), and by the pairs of running waves 

J±(t±s/v) = J 1 (s,t)±GV 1 (s,t), V±{t±s/v) = V 2 (s,t)±RJ 2 (s,t) 

of current and voltage in the first and second lines respectively as the solutions to the 
telegraph wave equations. The boundary conditions V\ (0, t) = V(t), J2(0, t) = ,J(t) 
induce on the open oscillator the following pair of Langevin equations 

CV{t) + GV(t) - J{t) = J+ (i) , LJ(t) + RJ{t) + V{t) = V+ (t) , 

J_(t) = J + (t) - 2GV(t) , V-(t) = V+(t) - 2RJ(t) , 

where the second pair of equations, corresponding to these boundary conditions, 
determines the output waves running from the oscillator to the receiver. 

Assuming, for simplicity, that the following gage invariance condition R/L = 
G/C = 7 is fulfilled, the above system of equations can be written in the one- 
dimension complexified form 

x(t) + (iü+^-f\x(t) =x+(t), x-(t) =x+(t) -7^(í) 

in terms of the complex amplitudes 

x(t) = ^(V(7F(t)+h/Zj(í)) , 

I±(i) " 7ï(75 J±{t) + 7z v±lt) y 

Now we move from the Langevin towards a stochastic description of the cor- 
responding quantum system. In the classical case, assuming that the input noise 
x+ (í), radiated by the receiver at the time t — s/v, is the equilibrium noise of the 
temperature T, we can treat it as a complex circular-invariant white noise with the 
intensity given by the Nyquist formula a = kT (where k is the Boltzman constant). 
In the quantum case complex amplitude x(t) should be replaced by the correspond- 
ing operator x(t) in the Heisenberg picture, satisfying the canonical commutation 
relations 

(2.1) [x(t), x(t)*} = x(t)x(t)* - x(t)*x(t) = hl , 

where x(t)* is the Hermitian adjoint operator and 1 is the unit operator in the 
Hilbert space of quantum-mechanical states of this system (h is the Planck constant, 
x(t)/y/h = à{t) is the annihilation operator for oscillator quanta at the time t). 
By virtue of linearity the corresponding Heisenberg equation for x(t) should have 
the same form as for the classical complex amplitude x(t), where the amplitude 
of the propagating wave x+ have to be replaced by the operator x+(t) with the 
commutation relations ensuring preservation of the canonical commutation relations 
(12. 1|) . This can be achieved if we take 

(2.2) [x±(t),x±(t')]=Q, [x±(t),x±(t')*]=' Y h6(t-t'), 



4 



V. P. BELAVKIN 



where we have taken into account that the commutators for the output wave X- (t) 
should coincide with the commutators of x+(t) due to its linear relation with 
x+(t), as it is also assumed that commutates both with x(r), x_(r) and with 

x(r)*, i-(r)* for r < t. As shown in [d], the above equations and commutation 
relations can be obtained by a canonical quantization of any open oscillator in the 
"rotating wave" representation and the narrow-band approximation. Introducing 
the notations 

(2.3) v(t) =x + (t) = -v(t), y(t)=x-(t-s/v) = -v(t) 

for the forward and backward waves at the input and output of the transmission 
line of length s respectively, let us write these equations in the following Standard 
form 

(2.4) f(f) + ax(t) = v(t) , y{t + s/v) = 1 x{t) + v(t), 
dt 

where a — iQ + 7/2. The first equation, rewritten in terms of quantum stochastic 
differentials as 

dx(t) + ax{t)àt = dv t , 
where Vt = f v(t)dt, can be easily integrated. The solution 

x(t) = e- Qt í(0) + f e a(r_í) dfi r 
Jo 

does not depend on the output wave y(t') with t' < t + s/v and it commutes with 
y(t') as well as with y(t')* for all such t'. This commutat i vity, reflecting the ideal 
conjugacy between the output measurement device and the transmission lines, will 
play an important role for the causal prediction of x(t), called the nondemolition 
condition. This condition, together with Markovianity of x(t) corresponding to 
the assumption that v (t) is a quantum white noise, which is fulfillcd under the 
narrow-band approximation for the equilibrium state of any measurement device 
sending the quantum thermal noise wave v(t) into the line towards the oscillator, 
significantly simplifies the optimization problem of the continuous measurement in 
the quantum system under the consideration. 

Let us now consider the problem of nondemolition observation of the amplitudc 
operator x(t) of the quantum oscillator at the output of the quantum noisy channcl 
described as above. One would like to obtain an optimal estimate x(t) of this 
amplitude as a classical stochastic complex-valued process, being observed in a 
real time on the output of the transmission line. The mean square optimization 
problem defining at each t such directly measurable nonanticipating estimate x(t) 
based on the results of the previous, in general indirect quantum nondemolition 
measurements, is the minimization problem of the mean quadratic error 

(2.5) ((í(t)-a:(í))*(*(í)-s(t))). 

The continuous in time quantum nondemolition measurement is defined by any 
measurable non-anticipating transformation of the received quantum process y (t) 
into the classical complex-valued one x (í), which should be in general randomized 
in order to include also indirect nondemolition measurements of the noncommuting 
y(t)* and y(t). Identifying the classical process {x (t)} with a family of normal 
compatible operators in an extended Hilbert space, one can describe it at each 
t > by a quantum stochastic functional x(t, •) of the past output operators 
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y(t')*, y{t') 1 if < t and also of some additional independent quantum noises ran- 
domizing if necessary this estimate, with the operator value x(t) which commutes 
with all operators x(t')* , x(t') even if t' > t. Note that such estimator x(t) will 
not necessary commute with all past estimated operators x(t' ,x(t')* , although due 
to the nondemolition property it will commute with them and all their future if 
t' > t — s/v. Such quantum directly observable process x(t), which is the closest 
in the mean square sense J22J) to the non-observed quantum stochastic process 
x(t), can be considered as a dynamical realization of the optimal temporal quasi- 
measurement (or estimation) of the noncommuting x(t' , x(t')* , based on the indirect 
observation of the given output process y (t) . 

In order to demonstrate such optimal temporal quasi-measurement let us con- 
sider the problem of quantum filtering in the Gaussian case when the initial quan- 
tum oscillator state is a circular-Gaussian with zero mathematical expectation 
(x(0)) = and average number of quanta E > 0: 

(2.6) (x(0) 2 ) = , (x(0)*x(0)) = hY, , 

and the quantum noise v(t) is also circular-Gaussian white noise: (v(t)) = 0, 

(2.7) («(í)ü(í')) = , (v(t)*v(t')) = h vl 5{t - t') , 

where v = (exp{?r//À;T} — 1) _1 is the average number of quanta in an equilib- 
rium state of the receiver with the temperature T. As will be shown in this paper, 
the optimal estimate x(t) of the operator x(t) mininhzing the quadratic criterion 
(|2.5(l is given by a nonorthogonal Hermitian-positive Gaussian operator- valued mea- 
sure which describes the coherent quasimeasurement of an integral nonanticipating 
transformation of the output process y. It can be realized by the direct measure- 
ment of the classical complex random process x(t) described by the linear Langevin 
equation 

(2.8) — x(t) + ax{t) = n{t) (y(t) - 7 x(í)) 

with x (0) = (x(Q)) — 0. The input process y(t) in this equation is a classical 
complex (i.e. commutative normal) process given by the additive randomizing 
transformation y (t) = y(t) + w(t) of the received quantum process y(t), where w(t) 
is the complex amplitude of an additional independent quantum noise with the 
opposite commutators 

[w(t), w(t')} = , [w(t), w(t')*} = -hjó{t - t'), 

zero expectations (w(t)} = and the minimal vacuum state correlations 

(2.9) (w(t)w(t')} =0, (w(t)*w(t')) = hjS(t-t') . 

Since y (i) is equivalent to the classical <5-correlated complex process of the intensity 
cr = H (y + 1) 7, this Langevin equation can be identified with the classical Kalman- 
Bucy filter which is usually written in the stochastic differential form as 

dx(t) + ax(t)dt = n(t)(dy t - jx(t)dtj , 

in terms of the continuous diffusive input process y t = J * y{t)àt. In this equation 
k(í) = (E(í) — v) /(l + v) , E(í) is the solution of the Riccati equation 

(2.10) Jl s(í) = ^< _ s(í))(1 + E(í)) ( 
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with S(0) = So and the effective classical white noise intensity a = Hj(l — 
exp{— /ry/fcï 1 }) -1 Thus the randomized nondcmolition quasimeasurement of x (t) 
is realized by the direct observation of x (t) in the result y (t) of indirect temporal 
observation of y (í) by the heterodyning processing y (t) i— * y{t) + w(t) Q, whcrc 
w(t) plays the role of the reference quantum wave. The mean quadratic error (12.51) 
for optimal filtration in this case has the minimal value ftE(í), which is not zero 
even under zero temperature T = corresponding to a = fvy, in the sharp contrast 
with the classical case when a — at T = 0. This corresponds to adding into the 
measurement channel a vacuum quantum noise which makes possible the hetero- 
dyne indirect measurement of y (t) statistically equivalent to a classical filtering of 
the complex white noise of the minimal intensity tvy > 0. Note that this conclusion 
remains also vàlid for the quasi-classical oscillator corresponding to Eo 3> 1, being 
indirectly observed in a quantum wave line. However in the infinite temperature 
case T — > oo when a 3> 1, the quantum consideration does not produce a substan- 
tial increase of the optimal filtring mean quadratic error compared with the one 
given by the classical Kalman-Bucy filter. 

3. Causal quantum measurement processes and filters 

Since the sequences x(ti), . . . ,x(t2) of the discrete-time measurements for arbi- 
trary {t\ < . . . < t n } completely determine a continuous function x(t) under the 
temporal observation, the quantum measurement process can be statistically de- 
scribed by the multitimc probability distributions P(da;(íi), . . . , àx(t n )), n = 1, 2, . . . 
which define the probability measure P (d{a;(í)}) as a projective limit on the space 
of the trajectòries {x(t)}. The statistical structure of quantum theory requires the 
probability measure to be a linear form with respect to the density operator g of a 
quantum state, represented as 

(3.1) P(d{x(í)})=Trèn(d{x(í)}) 

Here Tr means the Hilbert space trace, II (d{x(t)}) > is a Hermitian-positive 
measure on the space of the observable trajectòries {x(t)} with vàlues in an output 
operator àlgebra B. We call such IT the operator- valued probability measure (OPM) 
since it must be normalized to the identity operator J II (d{x(t)}) = 1 due to 
the normalization of P. This measure defincs, in particular, the multitime OPM 
ü(da;(íi) . . . dx(t n )) for any finite n = 1, 2, . . . which induce the probabilities of the 
time-discrete temporal observations as the linear functions of g according to l|3.1|l . 

Evcry decomposition of the identity 1 = J II(d{a;(í)}) into the positive operators 
II(d{x(í)}) corresponds, according to Naymark's theorem jS], to a measurement of 
some compatible set of operators, which are observed in an "extended quantum 
system" . But only those measures II correspond to the measurement processes 
physically realizable in real time, which satisfy the causality condition II* (dec*) € B 1 , 
where II* is the measure II on the space of "reduced" realizations x* = {x(s)} s <t, 
and {B f }, t > is a nondecreasing family of operator algebras B s Ç B*, defining the 
nondcmolition observables up to the time instants t. Typically B l is generated by 
the subfamily b f = {b(s) : s < t} of a given operator family {b(t)}, t > 0, describing 
an output quantum stochastic process b(t) as in the example of first Section where 

In order to have predictable behavior of a quantum object described by the 
Heisenberg operators {x (í)} generating an àlgebra A under the indirect temporal 
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measurements in the àlgebra B, the subalgebras B l must satisfy only the commuta- 
tivity condition B l C A' t for each t with respect to the algebras At generated by the 
present and future Heisenberg operators Xt = {x(s) : s > t} of the object such that 
there exist conditional expectation on At V B l with respect to £>' for each t. The 
measurement process should be also selfpredictable, for which we shall assume the 
existence of the conditional OPM n 4 .^, dx' ) on x\ = {^(í')}t'e(s,í]j describing by 
the formula 1)3. 1|) the non-anticipating processes of observation from any moment 
s < t up to t with known results x s = {x(t)}t< s of the preceding measurements. 
They should satisfy the compatibility condition 

(3.2) n'^Ad^n^Nd^) = UlKx'^dxll) 

for all t < tx < t 2 such that the Hermitian-positive operators II* (x s , dxl)* = 
n*(a; s , dec' ) commute II S (dx s ) € B s , having the vàlues in the relative commutants 
Bl of the output subalgebras B s with respect to B l . Such a compatible family 
{II^} s<t of the conditional OPM II^(a; s ,da;') £ B\, normalized to the identity oper- 
ator 1, will be called non-anticipating quantum filter. The measures II* , which are 
independent of x s , correspond to the filters without memory. The Markovian hlters 
described by the conditional measures II' , which depend only on the last preceding 
value x(s), are the simplest filters with the memory. In addition, the transition 
operator-valued measure U t s (x(s),dx(t)) defines, according to 1)3. 2JI . the multitime 
conditional measures 

Ií(x(t ),dx(t 1 ) ■ ■ ■ dx(t n ^i)dx(t n )) 
= Ull {x(t ),dx(h)) . . . II*«_ 1 (a;(t»_i), dar(í„)), 
satisfying the operator-valued Smolukhowsky equation 

(3.3) í n^( a; (ío),dx(í 1 ))n^(a ; (í 1 ),dx(í 2 )) = n^(x(ío),dx(í 2 )). 
Jx tl 

In the case of a linear measurable space X, for example a complex n-dimensional 
space C n , it is convenient to describe the Markovian filters by transitional operator- 
valued characteristic functions (OCF) defined as operator-valued Fourier integrals 

(3.4) $ï(x(a), u(t)) = J e i(u( ' )+x(t)+x(t)+u( ' )) n^(x(s), dx(í)), 

where u, x arc rt-dimensional complex columns, and x + , u + are their conjugate 
rows. The normalization condition for II' gives $ t J (x(s),0) = 1, and the condi- 
tion of Hermitian positivity II' > implies Hermitian positive-definiteness of the 
operator-matrix [<ï>(x(s), — u/)] > 0, where {uí,í = 1,2, . . .} is any finite col- 
lection of vectors £ C™. The operators $' commute on the non overlapping 
time intervals, and the Markovian condition 1)3. 2JI can be written in the following 
diffcrential operator form 

(3-5) ^-i^-y,u(í 2 )) ^(^Mh^Ut^o = ^(x(t ),u(í 2 )). 

Under certain regularity condition it can be proved that the opposite is also true: for 
every family of regular operator-valued functions {<£' , s < t}, which are continuous 
at u(t) = 0, satisfying the above normalization, Hermitian positive-definiteness and 
Markovianity condition 1)3.5)1 . the exists unique Markovian filter having the OCF 
$' 
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As ari example of a Markovian filter we can take a linear coherent filter which is 
defined by OCF function 

(3.6) $*(x(s), u(í)) = e iu(i)+í(t) e ii(í)+u(i) , 

where x(í) is vector-column composed of the operators x\(t), . . . , x n (t) and x(í) + 
is vector-row composed of the adjoint operators x\(t)* , . . . , x n (t)* . The operators 
x(í) entering into this anti-normal ordered expression for $1 are assumed to satisfy 
the following linear quantum stochastic differential equation 

(3.7) -^x(í) + B(t)x(t) = K(t)b(t), x(s) = x(s) 
dt 

with b(í) defined as m -dimensional column composed of the annihilation operators 
&i(í), . . . ,b m (t) of the quantum output process satisfying the following canonical 
commutation relations 

k(t),h(t')] = 0, [S fc , (í), 6,* (í')] = S H S(t - t 1 ), 

Here B(t) — [/3jw(í)], K{t) — [kjw(í)] are complex matrices of the size n x n and 
n x to respectively, which may continuously depend on time. It is not hard to 
verify by integrating the equation 1)3.7(1 that the characteristic operator-function 
(|3.6(l satisfies all the above mentioned properties, including l|3.5|l . In fact it can 
be easily seen that, since the operators x(r) are linear transformations from the 
annihilation ones b(t), that í I3.6fl is the characteristic operator-function of the 
operator-valued Dirac 5-measure 

n* (x( s ), dx(í)) = N-{s (*(t), dx(t))} , 

which corresponds to the anti-normal ordering "AA_", when the creation operators 
b*(í) act to the left before b(t). It is well known that such measure is generated by 
the coherent projectors and the Lebesgue measure on the complex space C™. For 
example, for the case of a fixed x(0) = the measure II(í, dx) = Üq(0, dx) has the 
form 

(3.8) II(í, dx) = |í, x) (í, x|dA(í, x) , 
where |í, x) are the normalized right eigen vectors 

(3.9) x(í)|í,x) = x|í,x), xeC" 

for the operators x(í) defined by the equation ( 13. 7|) with the initial condition 
x(0) = 0. The element dÀ of the volumc in C™, normalizing the expression l|3.8|l . 
is given by 

dA(í,x) = — f~^77aï n" =1 d Re x»d Re x t , 

where C(t) is the Hermitian matrix of commutators [3è k (t) , xi(t)*\ = (C(t))ki sat- 
isfying the equation 

(t) + B{t)C (t) + C{t)B (í)+ = K(t)K (t) + 
with C(0) = corresponding to x (0) = 0. 
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A physical realization of the above linear coherent quantum Markovian fïlter as 
optimal quasi- measurement of the operators {x(i)} is given by a precise direct mea- 
surement of the Markovian process x(í) described by the linear classical stochastic 
differential equation 

(3.10) ^x(í) + B(í)x(í) = K(t)b(t) , b(í) = b(í) + c(í) • 

Here c(í) is independent of b(í) quantum "vacuum" vector noise with the opposite 
commutators 

[c k (t),c l (t')]=0, [c k (t) 1 c l (t')*] = -8 kl 5{t-t l ) 1 

which make the components of the sum b(í) commuting and also commuting with 
the components of b(t)*, and (c(í)) = 0, 

(c k (t)ci(t')) =0, (c k (t)*ci(t')) =S k iS(t-t'), 

and so (c k (t)ci(t')*) = 0. Therefore the momenta of the process x(í) coincide with 
the anti-normal momenta for x(í), in particular (x(í)) = (x(í)), 

{x k (t)x t {t)) = (x k (t)xi(t)) , (jc fc (í)x,(t)*) - (x k (t) Xl (ty) . 
4. Quantum Gaussian Diffusion and Linear Filtering 

Let x (í) be an n-dimensional quantum diffusive process with zero initial expec- 
tation (x(0)) = and the initial corrclations 

(4.1) (x k (0)xi(0)) = , (x k (0)*xi(0)) = {R )i k , 

where i?o is a given positive-definite matrix. The process is defined by the stochastic 
equation 

(4.2) Ax(í)+A(í)x(í) = J(í)à(í), 

where A(t), J(t) are matrices of size nxn and nxm respectively, and à(t) is whitc 
quantum noise satisfying canonical commutation relations 

(4.3) [à fe (í),a;(í')] = 0, [à k ,(t),àt(t')]=5 kl 6(t-t') 

with zero expectations (à(í)) = and nonzero normal correlation matrix Q(t): 

(4.4) {à k {t)ài{t')) = , {à k (t)*ài(t')) = (Q(t) lk S(t - t'). 

We consider also an output system described at the output of a quantum linear 
noisy channel 

(4.5) 6(í) = F{t)±{t) + a(t) 

by to x n matrix F(t) and by m-dimensional quantum white noise à(í) with canon- 
ical commutators 

(4.6) [à k (t),à l (t')} = 0, [à k (t),à l (t')*]=S kl S(t-t'), 

zero expectations (à(í)) = and a nonzero normal correlation matrix N(t): 

(4.7) (à k (t)ai(t')) = , (a fc (í)*ai(í / )) - (N(t)) lk S(t - t 1 ). 

The pair (à(i),à(t)) is assumed to be independent of x(0) at least in the wide 
sense, having zero correlations with x(0)*, but it can have nonzero matrix T{t) of 
the mutual normal correlations in 

(4.8) (à k (t)ài(t')) = , (a fe (í)*a;(í')) = (T(t)) kl S(t - t'). 
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This is because we cannot assume them independent but satisfying the commutation 
relations 

(4.9) [fi fc (i), fi,(t')] = , [à k (t), ài{t')*\ = {D(t)) kl 5(t - t'), 

which are necessary for the nondemolition condition of mutual commutativity of 
all components b(í') with x(í)* for any t > t' . If x(í) satisiies the commutation 
relations 

(4.10) [x k (t),x l {t)]=0, [xk{t),xi{t)*] = (C{t)) H , 

the matrix D is defined by the nondemolition condition as a solution to the equation 
JD + CF + = at cach t (As it oceurred in the example considered in the Section 
121) Usually the commutation relations C (t) of the object are preserved, C (t) = Co, 
i.e. ACq + CqA + — JJ + . Note, however, that the nonzero commutation relations 
in (|4.3|) are not necessary for this preservation and for the nondemolition property 

[x k (t),y l (t')}=0, [x k (t)Mt'T] = Ví>t' 
if the process x(í) has the degenerate matrix Co, e.g. if it is classical: Co = 0, in 
which case also D = for the nondegenerate J. 

The problem of optimal nonstationary quantum filtering consist of finding a 
nonaticipating physically realizable measurement process, denoted by x(í), at the 
output of the quantum channel, minimizing at each time t the mean of the quadratic 
error 

n 

(4.11) |x(í) - x(í)| 2 := " x^fim " 

2=1 

Theorem 1. Let the initial quantum vector x(0) be Gaussian with zero mathe- 
matical expectation and the correlations \4-. the pair (à(í),a(í)) be also Gaussian 
with zero expectations and with correlations i4-7[ ), i4-8j ), ond a(í) satisfy the 

canònica! commutation relations \4-ty - Then the optimal by the criterion \4- 1 1\ ) 
quantum filter for the quantum Markov Gaussian process x(í) , defined by the linear 
quantum stochastic equation \4- l· ^il , is the Markovian coherent filter 1^3. b}) . described 
by the linear quantum stochastic equation \3. 7| ), where x(0) = and 

(4.12) B = A + KF, K = (PF+ + JT+)(N + I)- 1 . 

Here I is the identity m X m-matrix, and P = P(t) is the a posteriori correlation 
matrix, which is given by the solution of the Riccati equation 

(4.13) -^P + AP + PA + + (PF+ + JT + ){N + /) _1 (FP + TJ) = JQJ+ 

with the initial condition -P(O) = Rq. The minimal expectation value of the qua- 
dratic error \4-H\j at the time t is defined by the trace of P(t). (The dependence 
of all the matrices in \4-l!3j) , \4- 13{ ) on t is omitted for brevity.) 

Proof. Given the fact that the criterion for the point filtering l|4.11|l depends on 
the value x(í) of the observation {x(í)} only at the last instant í, the average value 
of the quadratic error is defined only by the single-time operator-valued measure 
II(í, dx(í)). Therefore one can find the solution II (í, dx(í)) of the statistical prob- 
lem ^Ü] of the optimal observation independently for each t. If it proves to be 
a single-time operator-valued measure for some physically realizable measurement 
process II (d{x(í)}) in the sense of the Section 03 then the problem of optimal 
filtering will be solved. 
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Using the stochastic calculus for normal ordered quantum Gaussian variables, 
which was developed in one can obtain the following representation for the 
average of the quadratic error (|4.11|l as a function of the operator-valued measure 
n(i, dx): 

n „ 

(4.14) (|x(í) -x(í)| 2 ) =J2(P(t)h + I Triï(í,x)n(í,dx). 

2—1 

Here P(t) is a matrix satisfying the Riccati equation (|4.13fl . R(t, x) is the operator- 
valued function 

n 

(4.15) R(t,x)=^(x i (t)-x i yè(t)(xi(t)-Xi), 

i=l 

of the operators x(í) satisfying the equation l|3.7fl defined by the matrices (|4.12f) 
with the initial condition x(0) = 0, and g(t) is a Gaussian density operator for the 
output process x(í). This g(t) can be represented in the form of a normally ordered 
expression 



g{t) = det(G(R - P)- 1 )^ 



k.l 



where G(t), R(t) are the solutions of the linear equations 

^-G + BG + GB+ = KK+, ^-R + AR + RA+ = JQJ+ 
dt dt 

with the initial conditions G(0) = and R(0) = Ro respectively. 

The second element of the right-hand side in (|4.14() is not negative because of the 

Hermitian positivity of the operator (|4.15|) , and it is equal to zero on the coherent 

measure (I3.8|) defined by equation i|3.9l ) for x = 0. Hence, it is clear that the 

coherent measurement, described by this measure, is optimal for all t. Therefore, 

the coherent observation described by the OCF $ (t, u) = $q (0, u) in with 

the initial condition x(0) = is the optimal measurement process. The Theorem 

is proved. 

Note that the a posteriori state for the quantum diffusion x(í) is described by 
the Gaussian stochastic density operator 

g(t) = det(G(P + G)- 1 )^ [exp {- " a*(t))*(P + G)w (*J " *«(*))} 

given by the solution x (í) of the classical stochastic diffusion equation l|3.10|l with 
x(0) = 0. 

As example, let us consider the n-dimensional quantum oscillator (|4.2|l that has 
à(t) = à(t), JJ + — h(A + A + ), such that the commutators 

(4.16) [xk,xi] = 0, [xk,x*] = hóki 

are preserved. The matching m-dimensional quantum communication line is de- 
scribed by the output wave 

(4.17) y(t) = G(t)±(t) - J(í)à(í). 

This wave commutes for G = A+A + with x(í')*, t' > t. In the former dimensionless 
units this line is described by the wave 1)4.5(1 . where F = —J + /h, so that y(í) = 
hF + h(t). Let the correlation matrix N of the quantum channel noise à(í) be 
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proportional to the identity matrix: N = vi. Then the equation (|4.13|) for the 
dimcnsionless a posteriori correlation matrix S = P/h will have the form 

(4.18) S = i[S, H] - ((S + I)G(S - vi) + (S- vI)G{S + /)) , 

where H = (A - A+)/2i, and K(t) = (ui - S)J/(1 + v). In particular, in the 
one-dimensional case m = n = 1, denoting S — S, G = 7, we obtain l|2.1ü|l . It 
is interesting to note that the stationary solution S = vl of the equation (|4.18|) 
corresponds to a nonstochastic optimal filtering which does not require at all any 
observation since in this case the amplification coefficient K(t) is zero. The initial 
condition Sq = vi for the stationary solution S (t) — So of the Riccati equation 
(|4.18|l in the case H = 7/ can be interpreted as a condition of thermal equilibrium 
for the initial statistical state of the oscillator having H = íl in the one dimensional 
case, and the the quantum channel noise equilibrium state as it corresponds to 
equality of the temperature for the equilibrium state in the communication line 
and of the temperature T = M7/Mn (l + for the oscillator initial Gaussian 

with E — 7- 
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